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The dynamical phase diagrams for vortices in clean films, driven by an uniform force, and inter- 
acting with periodic pinning resulting from a columnar-defect lattice are investigated by numerical 
simulations of a London model, and other considerations. Dynamical phases are identified accord- 
ing to the spatial symmetry, and dynamical phase diagrams are studied as a function of the driving 
force magnitude and direction, and on the temperature, for typical vortex densities. The theoretical 
limit of infinite drive- force magnitudes is investigated first. In this limit vortex motion averages the 
periodic pinning potential in the direction of motion, and the dynamical phases reduce to thermal 
equilibrium ones for vortices interacting with the averaged pinning potential. For most directions 
of motion this potential is essentially constant. For a few particular ones it is a one-dimensional 
washboard, periodic in the direction transverse to motion. The dynamical phase diagrams in this 
limit are obtained as a function of the direction of motion and temperature by equilibrium Monte 
Carlo simulations. They contain moving vortex lattices, commensurate or incommensurate with the 
washboards periodicities, at low temperatures, moving smectics and moving liquids at higher ones. 
Finite-drive dynamical phase diagrams are then obtained by numerical solution of Langevin's equa- 
tions. It is found that each dynamical phase originates from an infinite- drive limit phase with the 
same spatial symmetry that evolves continuously into finite-drive regions of the dynamical phase dia- 
gram. Dynamical transitions between moving commensurate or incommensurate lattices and moving 
liquids and between moving commensurate lattices and moving smectics are observed. Transverse 
pinning of vortices at low temperatures is found in regions of the dynamical phase diagrams where 
the dynamical phases are moving commensurate lattices or moving smectics. 
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I. INTRODUCTION 



A problem of current interest is the study of moving 
vortices interacting with arrays of pinning centers. Ran- 
dom arrays have received the most attention. However, 
periodic pinning arrays are also of interest. One rea- 
son is that they provide examples of dynamical phases 
and phase transitions, a subject not yet fully under- 
stood, where theoretical predictions can be tested in su- 
perconducting films with artificial defect latticestl, and in 
Josephson junction arrays (JJA)B. 

In the case of moving vortices interacting with random 
pinning arrays, the dynamics is simpler when the driving 
force is large, in which case the velocity of the center of 
mass (CM) of the vortex array is also large. This, situ- 
ation was first considered by Schmid and HaugerB who 
suggested that fast moving vortices average out the pin- 
ning potential and, consequently, order in a [triangular 
lattice at low temneratures. In experiments^ and nu- 
merical simulationsQ nearly triangular lattices are indeed 
observed at large CM velocity. Koshlev and VinokurB 
analyzed the effects of random pinning at large CM ve- 
locities on vortex displacement from the lattice equilib- 
rium positions. They argued that it is equivalent to a 



"shaking" temperature, inversely proportional to the CM 
velocity, that adds to the thermodynamic temperature, 
and predicted that dynamical melting of the moving vor- 
tex lattice takes place when the CM velocity is such that 
the combined temperature equals the equilibrium melt=. 
ing temperature. Recently, Gianmarchi and Le DoussalQ 
pointed out that this picture for the large CM velocities 
behavior breaks down, because the vortices order in a 
moving glass, rather than in a moving triangular lattice. 
The reason is that averaging of the pinning potential by 
the fast moving vortices is only partial. A static ran- 
dom pinning potential remains acting on the vortices in 
the frame of reference moving with the CM velocity (CM 
frame), and has non-trivial consequences on vortex or- 
der. These authors predict that the vortices move along 
static channels, on average parallel to the direction of 
drive, and are pinned on them with respect to a trans- 
verse force at zero temperature. This is referred to as 
transverse pinning. These predictions ape in agreement 
with experimental and numerical resultscl. 

In the case of vortices interacting with periodic pin- 
ning, the dynamics also simplifies in the limit of large 
driving forces. The question of what happens to the dy- 
namical phases in this limit was first examined in recent 
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publications by the authoru. These papers consider pin- 
ning by a square lattice of columnar defects (CD), and 
vortex motion in the [0,1] (or [1,0]) direction. In Refs. ||it 
is argued that moving vortices average the periodic pin- 
ning potential only in the direction of motion. It is also 
argued that the dynamical phases reduce to the equilib- 
rium ones for vortices interacting with the periodic pin- 
ning potential averaged in the direction of motion, that is 
[0,1], which is a washboard, periodic in the [1,0] direction. 
This was shown to have non-trivial consequences for the 
dynamical phase diagram. The present paper studies fur- 
ther the limit of large drives for the model introduced in 
Ref. ||, allowing for arbitrary directions of vortex motion, 
and covering a wide temperature range. The dynamical 
phases in this limit are found to play a central role in 
determining the whole dynamical phase diagram. 

Previous theoretical studies of vortices interacting with 
periodic pinning arrays are mostly numerical. Extensive 
studies of the zero-temperature vortex dynamics in sys- 
tems with CD-lattices, carried out by several workersO, 
find a rather complex behavior, with many dynami- 
cal phases and phase transitions. A simpler finite- 
temperature behaviop. is found in JJA and films. Mar- 
coni and DomimguezBI carried out numerical simulations 
of square JJA, with the driving force along the [1,0] (or 
[0,1]) direction . They find three dynamical phases: two 
moving vortex lattices and a moving liquid. In all three, 
the CM velocity is along the direction of drive. The lat- 
tices have the same spatial order, but differ from each 
other by the response to a force transverse to the direc- 
tion of motion. One lattice, named transversely pinned 
vortex lattice, is pinned (transverse pinning), and the 
other, named floating solid, is unpinned. Dynamical 
melting of the floating solid into the liquid, and dynami- 
cal depinning from the pinned lattice to the floating solid 
are reported. Further work£3 predicts transverse pinning 
only for driving forces along [1,0] ( or [0,1]), and verify it 
experimentally. 

WorkpQn films was carried-, out by Reichhardt and 
ZimanyiO and by the authora. Both investigate vor- 
tices interacting with a square CD-lattice driven along 
the [1,0], or [0,1], directions, and consider the following 
vortex densities: 1 vortex per CD in Ref. |ll|; 5/4 and 
2 vortices per CD in Ref. ||. They find that the direc- 
tions of vortex motion and drive coincide, and that the 
dynamical phase diagrams depend on the vortex density. 
For densities 1 and 5/4 vortices per CD the dynamical 
phase diagrams are similar, with three dynamical phases: 
a moving lattice, commensurate with the pinning poten- 
tial periodicity along the direction transverse to motion, 
a moving smectics and a moving liquid. Dynamical tran- 
sitions between the commensurate lattice and the smectic 
and between the smectic and the liquid are found. In the 
commensurate lattice with 1 vortex per CD tranverse pin- 
ning occursLLil. For 2 vortices per CD the only dynamical 
phases are an incommensurate lattice, essentially trian- 
gular, and 2. liquid, and dynamical melting between them 
takes placet! The present paper shows how these partic- 



ular film results fit into a dynamical phase diagram that 
depends on the magnitude and direction of the driving 
force, and on the temperature. 

The approach adopted here is similar to that used in 
Ref. ||. First the limit of high drives, where the dynami- 
cal phase diagrams are simpler, is studied in detail. Only 
then the paper proceeds to investigate the intermediate 
and low drives behavior. 

The paper starts by considering Langevin's stochastic 
dynamical equations for a generic model for vortices in- 
teracting with periodic pinning at temperature T, and 
its behavior in the limit where the driving force magni- 
tude approaches infinity. It is shown that in this limit 
these equations, written in the CM frame, describe vor- 
tices interacting between themselves and with a static 
effective pinning potential, equal to the average of the 
periodic pinning potential in the direction of motion. As 
a consequence, the dynamical phases reduce to the corre- 
sponding thermal equilibrium ones, which are hereafter 
referred to as infinite-drive phases. These phases are ex- 
act asymptotic limits of dynamical phases occuring in the 
dynamical phase diagram. The infinite-drive phases, be- 
ing equilibrium phases, are distinguished from each other 
by the spatial order. Thus, dynamical phases with dis- 
tinct spatial orders, one for each infinite-drive phase are 
expected to exist in finite-drive regions of the dynamical 
phase diagrams. The question is whether or not these dy- 
namical phases exhaust the dynamical phase diagrams. 
The answer is found be yes for the model studied in this 
paper, as well as for those investigated by Marconi and 
DominguezBI and by Reichhardt and Zimanyiliil. 

Next, the paper studies numerically the London model 
for vortices interacting with a square CD-lattice intro- 
duced in Ref. ||. Two vortex densities, typical of the 
region where there are more vortices than CD (B > B^), 
are considered: 2 vortices per CD (B = 2B$) and 5/4 
vortices per CD (B = 1.25£?0). The model dynami- 
cal phase diagrams for these two densities are obtained 
numerically. First, the model infinite-drive phases and 
phase diagrams are obtained as a function of the temper- 
ature and direction of drive by equilibrium Monte Carlo 
simulations of the London model for vortices interact- 
ing with the above mentioned effective pinning poten- 
tial. Then dynamical simulations of Langevin's equa- 
tions describing the model dynamics are carried out. 
Based on the results of these simulations, the dynami- 
cal phases are identified and characterized, and the dy- 
namical phase diagrams are obtained, as a function of the 
driving force magnitude and direction and of the temper- 
ature. The main findings are. i) All dynamical phases 
have a corresponding infinite-drive phase with the same 
spatial symmetry. The portion of the dynamical phase 
diagram occupied by each dynamical phase is found to 
extend continuously into large driving-force magnitude 
regions, where the dynamical phase is essentially identi- 
cal to the corresponding infinite-drive one. ii) Directions 
of drive and vortex motion are in general distinct, lead- 
ing to anisotropic V-I characteristics, iii) Existence of 
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two dynamical transitions. One, referred to as dynamical 
melting, between a moving incommensurate lattice and 
a moving liquid. Another, between a moving commensu- 
rate lattice or a moving smectic and a moving liquid, iv) 
Transverse pinning at low temperatures within regions of 
the dynamical phase diagram where the spatial symme- 
try is that of a moving commensurate lattice or a moving 
smectic, with vortex motion restricted to the [± 1,1] [0,1] 
and [1,0] CD-lattice directions. 

A brief summary of the main numerical results is given 
next. The CD-lattice and the definitions of the coordi- 
nate system, CM velocity and driving force orientations 
are shown in Fig. The model periodic pinning poten- 
tial is shown in Fig. ^.a. The effective pinning potentials 
resulting from it are washboards, if motion is in the [0,1], 
[1,0] and [±1,1] CD-lattice directions, and essentially a 
constant otherwise. The [0,1] and [-1,1] washboards are 
shown in Fig. and c. The [1,0] and [1,1] washboards 
follow from these by 90° rotation. The infinite-drive 
phases at low temperatures are found to be moving vor- 
tex lattices. For motion in the [0,1], [1,0] and [±1,1] 
directions these lattices can be commensurate or incom- 
mensurate with the washboard's periodicity, as shown in 
Fig. ||. For motion in other directions the moving lattice 
is essentially triangular. Dynamical phase diagrams for 
driving forces oriented in the [0,1] and [-1,1] directions, in 
which case the directions of motion and drive coincide, 
are shown in Fig. ||. The low-temperature dynamical 
phase diagrams as a function of the driving force mag- 
nitude and direction are shown in Fie;. | Within MIL 
regions in Fig. j| a moving commensurate lattice exists, 
with the vortices ordered in a nearly triangular lattice. 
Dynamical melting lines separate it from a moving liq- 
uid. In the regions denoted MCL MCL45 and MSM in 
Fig. |[ transverse pinning occurs, with vortex motion is 
restricted to the [-1,1] direction. The vortices order in 
moving commensurate lattices within MCL and MCL45, 
with spatial symmetries equal to the infinite-drive phases 
with the same labels in Fig. ||. Within MSM the spatial 
order is that of a moving smectic, similar to the MSM 
phase for B = 2B§ shown in Fig. [I[ 

This paper is organized as follows. The infinite-drive 
limit is considered in Sec. [|. The London model and 
the numerical methods used to study it are presented in 
Sec. QI. The simulation results are reported in Sec. IV. 
Finally, in Sec. [v| a discussion of the paper main results 
and comparison with earlier ones is presented, and the 
conclusions are stated. 



II. INFINITE-DRIVE LIMIT 

Here the limit of infinite driving force magnitudes is 
discussed for a generic model for vortices interacting with 
a CD-lattice. 

The motion of N v two-dimensional vortices is assumed 
to be governed by Langevin equations for massless par- 



ticles, which for the l-th vortex read; 



' dt 1 



j\v — cdl 
I 



where r] is the friction coefficient, 

N v 



- ViW-vfa-rj) 



(1) 



(2) 



is the force of interaction with other vortices, U v ~ v (r) 
being the vortex-vortex interaction potential in two- 
dimensions, 



■*v—cdl 



-ViU v - cdl (n) , 



(3) 



is the force of interaction with the CD-lattice, U v cdl (r), 
being the respective potential, given by 



U v ~ cd \v) = U v ~ cd (r - R) , 



(4) 



R 



where R denotes the CD-lattice positions and U v ~ cd (r) 
is the interaction potential between a vortex and a sin- 
gle CD, is the driving force, and Ti is the random 
force appropriate for temperature T. In terms of Fourier 



transforms U? cdl can be written as 



U 



v—cdl 



(r)=Vr d (Q)e'* r , 



Q 



(5) 



where Q denotes the CD-lattice reciprocal lattice vectors, 
and U v ~ cd (Q) is the Fourier transform of U v ~ cd (r). 

It is convenient to consider vortex motion in the frame 
moving with the CM. Let r[(t) = ri(t) — V cm t denote 
the Z-th vortex position in the CM frame, V cm being the 
CM velocity, to be defined shortly. In the CM frame the 
vortex-CD lattice interaction, denoted here by F^~ cdl (t), 
depends explicitly on time, namely 



F v-cdi^ = J2{-iQ)U v - cd (Q)e iQ ^e iQ - 
Q 



V crn £ 



(6) 



The equation of motion for the l-th vortex in the CM 
frame reads, 



] ~dt 



cdl 



(*) 



(7) 



The CM velocity is defined by the condition 
J2[ v {dr f l {t)/dt))t = , where () t denotes average over 
the random force distribution and over time. It follows 
from Eq. that 



V 3 = 1 



(8) 



In the limit Fd — > oo, it follows from Eqs. (g) and 
(§) that V cm = Fd/77, so that V crn — > 00 also. In this 
case, all Fourier components in F^ _cdZ (t), Eq. (||), for 
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which Q • V cm 7^ oscillate rfost, having a negligible ef- 
fect on the vortex trajectory^, and F^~ cdl (t) reduces to 
the static force obtained by summing the Fourier com- 
ponents in F[ v ~ cdl \t) with Q • V cm = 0, or Q 1 F d , 
namely 



■^v — cdl 
I 



(*) 



Q-LFd 



J (Q)e 



iQr'l 



(9) 



This force derives from the static effective pinning poten- 
tial 



K ff (r' ± ) 



(10) 



where r' ± is the CM frame coordinate perpendicular to 
the direction of F^, and a the orientation of F^. This po- 
tential is equal to the average of JJ v ~ cdl in the direction 
of motion (and drive). By definition, U^^(r f ± ) is one- 
dimensional and periodic in the direction perpendicular 
to Fd, if Fd is oriented along one of the square CD-lattice 
directions and a constant otherwise, since no Q is perpen- 
dicular to Fd- The equations of motion in the CM frame, 
Eqs. (0), reduce then to those describing vortices inter- 
acting between themselves and with U^/^{r f j_) at temper- 
ature T. Their long-time solutions are the corresponding 
equilibrium phases, called infinite- drive phases in this pa- 
per. 

The infinite- drive limit establishes the exact asymp- 
totic behavior Fd — > oo of the dynamical phase for given 
Fd orientation and T. For finite Fd, the dynamical phase 
is expected to be close to the infinite-drive one with the 
same spatial symmetry for sufficiently large Fd. The con- 
dition for this to happen is that the vortex displacements 
in the CM frame caused by the pinning force oscillations 
in time, Eq. (Q), during a time interval of the order of one- 
half period, are small compared with the vortex mean 
separation a v . This means that for every Q, such that 
Q • F d ± 0, 



Q | U^ v - cd \Q) 



<C a v 



(11) 



III. MODEL AND NUMERICS 

The model used in the numerical studies carried out 
here was introduced in Ref. |[ Vortices and CD are placed 
on a square lattice (space lattice), subjected to periodic 
boundary conditions with TV = 256 x 256 square prim- 
itive unit cells of dimensions d x d. The CD-lattice is 
square, oriented parallel to space lattice and commensu- 
rate with it. The CD-lattice, the coordinate system and 
the angles defining the orientations of F^ (a) and V cm 
(0) are shown in Fig. |l|. 

The vortex has a core of linear dimension d v = Ad 
(d v ~ 2£(0)), so that each vortex occupies 16 space- 
lattice sites. The vortex- vortex interaction potential 



is a screened Coulomb oneEl, periodic in the space 
lattice, and characterized by the energy scale J — 
((/>Q<i v /327r 3 A 2 ), where A is the penetration depth. That 
is, U v ~ v (r) is the lattice Fourier transform of 



U v ~ v (k) 



4tt 2 J- 



A- 



where n 2 = 4 sin 2 (k x d/2) + 4 sin 2 (k y d/2), k x and k y 
being the space lattice reciprocal lattice vectors compo- 
nents, k c = 2 sin {itd/2d v ) is the vortex core cutoff in 
k-space, and A is the screening length (A > A). The 
square CD-lattice has N c d = 8x8 sites, lattice constant 
a c d = 32d. The interaction potential between a vortex 
and a single CD is chosen with depth JJ v ~ cd (r = 0) = — J, 
range R c d = 12d and a spatial dependence that gives 
square equipotentials and a pinning force of constant 
modulus F p = J /R c d, as shown in Fig. |[a. The re- 
ciprocal CD-lattice vectors are 

Q = Qi(nix + n 2 y) , 

where Qi = 2tt / 'a c d, and ri\ ,ri2 = —4, —3,..., 3. It is 
found that the lattice Fourier transforms U v ~ cd (Q) are: 
U v ~ cd (0) = -0.19J, U v ~ cd {±Qi±) = U v ~ cd (±Q 1 y) = 
0.10J, U v - cd (±Q 1 [x±y}) = -0.06, and essentially neg- 
ligible otherwise. 

The model has the square lattice symmetry, so that 
the dynamical phase diagrams need only to be studied 
for Fd orientations 0° < a < 45°. 

The dynamical phases spatial order is obtained by 
calculating the time-averaged density-density correlation 
function, P(r), which is proportional to the probability 
that a vortex is found at r, given that there is one at 
r = 0, and its Fourier transform, the structure function, 
5'(k). The motion of vortices is characterized by calcu- 
lating the CM velocity, and the time-averaged velocity of 
each vortex. 

Two types of numerical studies of the model are car- 
ried out. Equilibrium Monte Carlo simulations to obtain 
the infinite-drive phases, described in detail in Sec. IV A . 
Numerical integration of Langevin's equations to obtain 
the dynamical phase diagram as a function of Fd, a and 
T for two typical vortex densities: two vortices per CD 
and five vortices per four CD. These corresponds to the 
magnetic inductions B = 2B$ and B = 1.25^, respec- 
tively, where B^ = </>o/a 2 d is the matching field, for which 
there is one vortex per CD. The numerical integration 
algorithm treats the vortex displacements as continuous 
during the run. In order to calculate the force acting 
on the vortices at each time step, the vortex positions 
are reduced to the nearest space lattice ones. The ac- 
curacy of this approximation increases with the number 
of space lattice sites occupied by each vortex, equal to 
16 here. The numerical integrations are carried out at 
constant 5, a and T, starting with a large Fd and the 
vortex positions initialized in a configuration appropriate 
for the infinite-drive phase. Subsequent runs, each with 
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a smaller Fd than the previous one, are carried out with 
the vortex positions initialized in the last configuration 
obtained in the previous run, until vortex motion can no 
longer be detected. In some cases, runs with Fd retracing 
the decreasing sequences are carried out in order to check 
for hysteresis. Typical runs consists of ~ 0.4 — 1 x 10 6 
integration steps at intervals 10 _2 r, where r = r]d 2 /J is 
the natural time scale. In order to average the thermal 
force, the data reported here is the average over 10 repe- 
titions of the same run, with different realizations of the 
random force. 



IV. RESULTS 



In this section the simulation results for the dynamical 
phase diagrams are reported. First in the infinite-drive 
limit, then for finite drives. In the figures presented here 
driving force magnitudes are measured relative to the sin- 
gle CD pinning force F p , Sec. Ill; temperatures relative to 
the infinite-drive movi ng inc ommensurate lattice melting 
temperature T m , Sec. IV A and center of mass velocity 
components relative to the components of V^ = Fd/rj. 



A. Infinite-Drive Limit 



The infinite-drive phase diagrams are obtained by equi- 
librium Monte Carlo simulations of the lattice London 
model with the pinning potential U%ff(r±). It is found 
that for F d along [0,1] (a = 0°) and [-1,1] (a = 45°) 
the U%ff(r±) are the washboards shown in Figs. [|b) and 
c). For Fd along other lattice directions, the U^/^(r±) 
are found to be very shallow washboards, because the 

uiv-cd)^ ( Eq> @) are very smaU ( Sec> and are 

considered as constant potentials in this paper. For ex- 
ample, the U%ff(r±) for along the [-1,2] direction (a = 
26.6°), shown in Fig. |^.d has a well depth more than one 
order of magnitude smaller than that for a = 0° and 45° 
shown in Figs. 0b) and c). For F^ oriented in non-lattice 
directions, U%f*(r±) is a constant. The corresponding 
low-T infinite-drive phases are as follows. For a = 0° and 
45°, they are vortex lattices (VL) commensurate or in- 
commensurate with the one-dimensional periodicity, de- 
pending on B. The VL for a = 0°, are obtained here as a 
function of B, for B > B^. Their density-density correla- 
tion functions, P(r), are shown in Fig. || The commen- 
surate lattices consist of identical vortex chains within 
each washboard channel, with neighbor chains displaced 
with respect to each other by half a chain period. There 
is a single chain for B^ < B < B\, and two chains for 
B\ < B < B2 . The values of B\ and B2 are found to be in 
the ranges B 1 < 1.125^ and 1.12550 < B 2 < 1.375P . 
A commensurate-incommensurate transition takes place 
for B = B c _i (1.25#0 < B c ^ < 1.375£ ). For B > B c ^ 
the VL is incommensurate, and nearly triangular. For 



B = 2P and B = 1.25P the P(r) shown in Fig. | cor- 
respond, respectively, to the vortex lattices labeled MIL 
and MCL0 in Fig. The effect of the one-dimensional 
periodic potential on the incommensurate lattice is to 
displace the vortices from the triangular lattice positions 
by a distance small compared to the lattice parameters. 
It is found that, as expected, this effect is greater for B 
in the vicinity of B c -i. This is evidenced by the behav- 
ior of P(r) shown in Fig. |. For B = 1.375P , 1.5P , 
and 1.7550, -P(r) has smeared spots centered in a nearly 
triangular grid, with smearing increasing with distance. 
For B > 2£>0, P(r) has sharp spots, indicating that the 
displacements are negligible, and the VL is sharply de- 
fined. 

The two B values studied in this paper, B = 2B§ 
and B = 1.255^, are typical of parameter regions far 
from commensurate-incommensurate transitions, both 
for a = 0° and a = 45°, where the VL are sharply de- 
fined. The infinite-drive phase diagrams in these cases 
are as follows. For a = 0° and 45° the VL are shown in 
Fig. j|. These are referred to in this paper as moving in- 
commensurate lattices (MIL) and moving commensurate 
lattices, with distinct vortex configurations labeled MCL, 
MCL0 and MCL45, as shown in Fig. |. For 0° < a < 45°, 
the VL are incommensurate and nearly triangular. It is 
found that the incommensurate lattices for the same P, 
but different a, cannot be distinguished from each other. 
Hereafter all incommensurate lattices are referred to as 
moving incommensurate lattices (MIL). 

The T-dependence of the infinite-drive phases (at con- 
stant a is found to be as follows: the moving incommen- 
surate lattices (MIL) melts into a moving vortex liquid 
(MLQ) at k B T m = 0.09J, for both B. For B = 2B^ 
the moving commensurate lattice (MCL) changes into a 
moving smectic (MSM) at T mc //T m = 1.4 and the moving 
smectic (MSM) changes into a moving liquid (MLQ) at 
Tmsm/Tm = 1.8. For B = 1.25^0 the moving commen- 
surate lattices, MCL0 and MCL45, change into a moving 
smectic (MSM) at T mcl0 /T m = 1.2 and T mcl45 /T m = 1.7, 
respectively. 



B. Finite Drives 

First the results for F^ oriented along the high sym- 
metry directions [0,1] ( a = 0°), and [-1,1] ( a = 45°) 
are presented. In these cases the directions of motion 
and drive coincide (0 = a) all the way down to the Fd 
values where V cm = 0. Next results for 0° < a < 45° are 
discussed. 



1. Drive along a = 0° and a = 45° 

Dynamical phase diagrams (Fd vs. T) are shown in 
Fig. ||. For B = 2P^, a = 0°, only two dynamical phases 
exist: a moving incommensurate lattice (MIL) and a 
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moving liquid (MLQ), separated by a dynamical melting 
line, as shown in Fig. [|. The dynamical phase diagrams 
for both B, a = 45°, and for B = 1.25^, a = 0°, are 
similar to one another, containing three phases: a mov- 
ing commensurate lattice (MCL for B = 2£?0, MCLO and 
MCL45 for B = 1.25^) a moving smectic (MSM) and 
a moving liquid. Only the dynamical phase diagrams 
for B 2£ , a = 45°, and B = 1.25£ , a = 0° are 
shown in Fig. |. That for B = 1.25£ , a = 45°, is 
similar. In the temperature range covered by Fig. [| the 
transitions between the moving smectics and moving liq- 
uids do not appear. The low T properties of the a = 
phase diagrams are studied in detail in Ref. ||. In Fig. 
U typical P(r) at large Fd/F p for the above described 
phases are shown. It is found that the dynamical phases 
at large Fd/F p essentially coincide with the correspond- 
ing infinite-drive ones. The dynamical phase diagrams in 
Fig. U show clearly that each dynamical phase originates 
from the infinite-drive one with the same spatial symme- 
try that evolves continuously into finite drive regions. 

For B — 2£?0 the vortices freeze at low T into the 
zero-drive lattice (ZDL in Fig. |l|) below the dotted lines 
in Fig. [I[ In this case the phase diagrams show little 
irreversibility. That is, if Fd is increased staring from 
the zero drive lattice essentially the same results are ob- 
tained. For B — 1.25^0 the vortices do not freeze at low 
T into the equilibrium zero drive lattice, but remain in a 
metastable state, similar to the moving smectic, down to 
very small Fd. The zero-drive lattice in this case is found 
to be a complex VL, commensurate with the full CD- 
lattice potential, but with many vortices per unit cell. 
However, if Fd is increased starting from the metastable 
state, little hysteresis is found. 

The transition lines show in Fig. ^ indicate the exis- 
tence of two regimes for the dynamical transitions. One, 
for large i 7 ^, where the transition temperature depends 
weakly on Fd, being close to its equilibrium values in 
the infinite-drive limit, showing that the transitions are 
driven by thermal fluctuations. The other, at small i 7 ^, 
where the transition force depends weakly on T, being 
close to the limit value of Fd as T — > 0, indicating that 
the transitions are driven by vortex vibrations, caused 
by motion in the periodic pinning potential. The lat- 
ter mechanism is similar to vortex 'shaking', proposed in 
Ref. [(] for random pinning, but a 'shaking' temperature 
cannot be defined. 

As discussed next, the commensurate VL ( MCL, 
MCLO and MCL45) show transverse pinning at low T. 
This may be expected because for a = 0° and a = 45°, 
V cm remains parallel to F^. In all these cases, according 
to Eq. (||), the effective pinning potential (U^^(r±)) acts 
on the moving vortices for all Fd , not only in the Fd — > oo 
limit. The commensurate phases are pinned by U%ff(r±) 
with respect to a small force along the direction defined 
by r_L. 



2. Drives along 0° < a < 45° 



The dynamical phases found in this range of driving- 
force orientations have the s ame spatia l symmetries as 
those described in Sees. [V A and IV B 1 and are referred 



to in what follows by the same nomenclature. 

The Fd vs. a dynamical phase diagrams at low-T are 
shown in Fig. ||. Both have two dynamical transition 
lines. One, referred to here as dynamical melting line, 
separating a moving incommensurate lattice (MIL) and 
a moving liquid (MLQ). Another separating a moving 
commensurate lattice (MCL for B = 2B§ and MCL45 for 
B = 1.25#0) or a moving smectic (MSM) for B = 2B^ 
and a moving liquid. Within the moving commensurate 
lattice and moving smectic regions transverse pinning oc- 
curs, with the vortices moving along the [-1,1] direction. 
The transition line in the B = 2B^ dynamical phase dia- 
gram from the moving commensurate lattice (MCL) and 
the moving smectic (MSM) occurs with the vortices mov- 
ing in the [-1,1] direction , and is essenti ally identical to 
that for a = 45° discussed in Sec. [V B 1 



For clarity of presentation, the dynamical phase dia- 
grams of Fig. [5] do not show the Fd vs. a line representing 
the cessation of vortex motion at small Fd. It is found 
that the behavio r is sim ilar to that for a = 0° and 45° 
discussed in Sec. |lVBll For B = 2B(f) the vortices freeze 
in the zero-drive lattice (ZDL) shown in Fig. [l] for all 

a. For B = 1.25£?0 freezing of the vortices occurs in a 
metastable state which change somewhat with a. 

The evidence leading to the construction of the dynam- 
ical phase diagrams in Fig. || is summarized in Figs. [?], 
||, H, and |HJ where the dependencies of the spatial order 
and V cm on i 7 ^, along lines of constant a in Figs. |[a and 

b, are shown. The evolution with Fd of 5(k) is shown 
in Figs. ^ and |s| along the a = 35° line for B = 25^, 
and along the a = 40° line for B — 1.25^. The changes 
occurring in S(k) shown in these figures when Fd crosses 
one of the transition lines are typical of those occuring 
across the same line at different a. In Figs. || and [l(] the 
direction of motion (0) vs. Fd curves are shown for var- 
ious a. Also shown in is the component of V cm perpen- 
dicular to the [-1,1] direction, denoted V±, as a function 
of Fd. This quantity is proportional to the voltage differ- 
ence between contacts placed perpendicular to the [-1,1] 
direction. In Figs. |^ and |l0[ V± is measured relative to 
the component of = Fd/r] perpendicular to the [-1,1] 
direction, denoted by Vd±. 

The detailed properties of the dynamical phase dia- 
grams shown in Fig. ^| are as follows, 
i) Dynamical Melting. The dynamical melting lines ex- 
tend from a = 0° to a = 45°. For B = 2B§ it touches the 
a = 0°-axis at the Fd value where the dynamical melting 
takes place for a = 0° and T/T m = 0.83 (Fig. |). It does 
not touch the a = 45° axis for both B and the a = 0° 
axis for B = 1.2550, because the dynamical phases in 
these axes are n ot m oving incommensurate lattices, as 
discussed in Sec. EVA . These result of (Fig. ||)show that, 
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for both £?, the moving incommensurate and commensu- 
rate lattices are separated by a moving liquid, at least up 
to the highest Fd studied here. However, in the infinite- 
drive phase diagram there is no moving liquid at these 
temperatures. It is unclear how the disappearance of the 
moving liquid as Fd — > oo takes place. 

The vs. Fd curves in Figs. ^ and [l(] show that in the 
moving incommensurate lattice (MIL) the directions of 
drive and vortex motion are approximately equal (0 ~ a), 
and that the melting transition is reflected in this curve 
by a change in slope. For B = 2B^ this change is sharp, 
with increasing rapidly towards = 45° as Fd de- 
creases. For B = 1.25£?0 the same is true, as long as 
a > 20°, while for smaller a the moving liquid slowly 
approaches = 0° as Fd decreases. 

It is found that the moving liquids are anisotropic, par- 
ticularly when is close to = 45°. This is evidenced 
by the S(k) panels in Fig. gfor F d /F p = 2.0, and in Fig. 
| for Fd/F p = 2.5 and F d /F p = 0.75, where a pair of 
peaks at k = ±Qi(x + y) occurs. These peaks result 
from a modulation in the vortex density along [1,1], with 
wavelength V^&cd, due to vortex motion being predomi- 
nantly along the [-1,1] direction (0 ~ 45°, see Figs. ^ and 
|l(J) , which favor concentration of vortices around lines in 
the [-1,1] direction connecting the CD-lattice potential 
minima (Fig. §.a). 

ii) TransverseP inning. The results shown in Figs. 
|9] and [l(] indicate that in the regions where the spatial 
order is that of a moving commensurate lattice, or a mov- 
ing smectic for B = 25^, transverse pinning occurs, and 
vortex motion is restricted to the [-1,1] direction. This is 
seen in Figs. ^| and [l(], where Vj_ = and = 45° in the 
Fd range of moving commensurate lattices or smectics, 
and the distribution of individual vortices direction of 
motion is sharply peaked around = 45° (Fig. ^|, inset). 
In the literature, tranverse pinning is usually discussed 
in terms of the component perpendicular to the di- 
rection of motion needed to depin the vortices. For the 
results shown in Fig. ^| this quantity, denoted here by 
F± c , can be obtained from the transition lines between 
the moving commensurate lattices, or smectics, and the 
moving liquid by Fj_ c = Fd sin (45° — a), where Fd and a 
are taken at these lines. These results also show that the 
moving commensurate and smectic regions evolve contin- 
uously from the co rresp onding infinite-drive phases. As 
pointed out in Sec. IV A , transverse pinning can be antic- 
ipated when the infinite-drive phase is a commensurate 
lattice. In the cases where this lattice is a single- chain 
one, a = 45° both B (Fig. ||), the depinning force can be 
estimated as being that for a single vortex, which is equal 
to the maximum force exerted by the a = 45° washboard 
( Fig. ^.c). The value of the maximum force is found to 
agree, within the simulation errors, with F± c obtained 
from the data shown in Fig. [5] at large Fd. In the case 
B = 1.25£?0, it is expected that a moving commensurate 
lattice region with spatial symmetry MCLO occurs near 
a = 0°. This region is not observed for T/T m = 0.89, 



used in Fig. g, but lower T results indicate that it exists, 
and is restricted to the immediate vicinity of a = 0°. 
This may result from a very small transverse depinning 
force for the MCLO moving commensurate lattice, due 
to its two-chain structure. This makes it difficult to ob- 
tain accurate results for this region, and it will not be 
discussed further in this paper. 

The temperature dependence of the dynamical phase 
diagrams is as follows. For temperatures lower than those 
of Fig. [5] the dynamical phase diagrams are similar. The 
dynamical melting lines are shifted to lower F^, similarly 
to what happens for B = 2B§ and a = 0° (Fig. ^ ). 
Transverse pinning in the moving commensurate lattices 
also occurs, as evidenced by the V± vs. Fd curves shown 
in Fig. |ll|. For higher temperatures several changes are 
observed. The dynamical melting lines also exist, shifted 
to larger F^ provided T does not exceed the infinite- drive 
melting temperature T m . When T > T m the moving in- 
commensurate lattices do not exist, only moving liquids. 

As well knownQ, transverse pinning in the moving com- 
mensurate lattices or smectics only exists, strictly speak- 
ing, at T = 0. At finite T thermal excitation of lattice 
defects lead to flux motion away from the [-1,1] direc- 
tion, giving rise to a finite V±. The regions of transverse 
pinning found in the present simulations occur because 
the thermally excited V± is smaller than the simulation 
resolution. In experiments, where the resolution is also 
finite, these regions may also be present. The T depen- 
dence of the V± vs Fd curves are shown in Fig. [ll|. For 
B = 2£?0 transverse pinning is found in the moving com- 
mensurate lattice and moving smectics at temperatures 
above T m , whereas for B = 1.25^ it disappears at lower 
temperatures. 



V. DISCUSSION 

The results described in Sec. |y| reveal the important 
role played by the infinite- drive phases. They show that 
all dynamical phases originate from infinite-drive phases 
with the same spatial symmetry that evolve continuously 
into regions of finite drives. _ 

The results obtained by Reichhardt and Zimanyitil for 
vortices in films interacting with a square CD-lattice, 
driven in the [1,0] direction, at B = B^, show a simi- 
lar behavior. The infinite- drive phase diagram was not 
considered in this paper, but can be identified with the 
reported dynamical phase diagrams at high drives. It 
consists, in the terminology adopted here, of the follow- 
ing phases. A moving commensurate lattice at low T, 
with single-chain structure similar to that shown in Fig. 
H for B = 1.1 251^ rotated by 90°, a moving smectic 
at intermediate T, and a moving liquid at high T. The 
dynamical phase diagram reported in Ref. 
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shows that 

these three phases exist in continuous regions of the driv- 
ing force vs. T phase diagrams that extend from high to 
low drives. This phase diagram is similar to those found 
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here for B = 1 .25^. a = 0° and 45°, and for B = 2B^ 



a = 45° ( Sec. |VB1|). 

For vortices in J J A, driven in the [1,0] direction, 
the dynamical phase diagram reported by Marconi and 
DominguezH also show dynamical phases existing in con- 
tinuous regions of the driving force vs. T plane that 
extend from high to low drives. However it is not pos- 
sible to identify the infinite-drive phases from the data 
reported in the paper. 

Some studies of vortices interacting with CD-lattices 
at T — 00 find dynamical phases and phase transitions 
similar to the ones reported in this paper and in Ref. pj] . 
Examples are tranverse pinning, moving commensurate 
and incommensurate lattices and dynamical transitions 
between them. However, these studies probe dynami- 
cal behavior different from the finite-temperature ones 
discussed here. At finite temperatures, the vortices re- 
lax after sufficient long times to a steady state, iden- 
tified as the dynamical phase, that is unique as far as 
the probability distribution is concerned!!^. Accordingly, 
dynamical phase changes caused by varying the driving 
force, the temperature, or both, are reversible. This is 
not necessarily the case at T = 0. However, in numeri- 
cal studies, irreversibility appears due to insufficient run 
times to reach the steady state, particularly when relax- 
ation times are very large. This occurs in several circum- 
stances, such as low temperatures, or vicinity of dynam- 
ical phase changes. In the simulations reported in this 
paper low temperatures are avoided. Some simulation 
runs where performed cycling the driving force magni- 
tude from a large value to a small one and back. Small 
hysteresis is found near dynamical phase transitions. 

In conclusion, the approach adopted in this paper to 
construct the dynamical phase diagrams, starting from 
the infinite- drive ones, provides a simple method, based 
on equilibrium statistical mechanics, to identify dynami- 
cal phases spatial order, and to predict dynamical phase 
transitions. It is expected that this method is applica- 
ble to a large class of periodic pinning potentials and 
vortex densities. The reasons are that the effective pin- 
ning potentials resulting from averaging physically rea- 
sonable periodic pinning potentials in the direction of 
motion are, as those discussed in Sec. IV A, essentially 



constant for most directions, and one-dimensional and 
periodic for a few particular ones. This predicts the ex- 
istence, for B > B^, of moving lattices, commensurate 
or incommensurate with the effective pinning potential 
periodicity, moving smectics and moving liquids, similar 
to the ones reported here, and of dynamical phase transi- 
tions between then. For B < B^ the moving lattices may 
be different, particularly at low vortex densities, where 
the commensurate infinite- drive phases are not expected 
to retain the simple chain structure found here for the 
B > B^ ones. However, the close relationship between 
the dynamical phase diagrams and the infinite-drive ones 
are still valid, as evidenced by the JJA results of Ref. ||. 
The dynamical phases obtained by this method do not, 
in general, exhaust the dynamical phase diagram. One 



known type of dynamical phase that has no correspond- 
ing infinite-drive phase is that in which some vortices are 
pinned by the periodic potential and others are moving. 
These are found at low drives in some T = simula- 
tions for B > BjB. No such phases are found in the 
present simulations nor in Refs. |],[)| The details of how 
the dynamical phases and dynamical phase transitions 
predicted by the method proposed here fit into the dy- 
namical phase diagram for each particular model and vor- 
tex density depends in a complicated way on the model 
parameters, and has to be determined in each case. 
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FIG. 3. Vortex positions in low-T infinite- drive phases. 
Dashed gray lines indicate minima of U%/^ shown in Figs. 2. b 
and 2.c. Nomenclature: moving commensurate lattices: a) 
MCLO, b) MCL45 and d) MCL. Moving incommensurate lat- 
tices: c) MIL 



FIG. 1. Columnar defect-lattice and definitions of coordi- 
nate system, CM velocity (V C m) and driving force ( F^) ori- 
entations. Inset: low-T equilibrium vortex state at zero-drive 
for B = 2B^ (ZDL). 
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FIG. 2. a) Columnar defect-lattice pinning potential, b-d) 
Effective pinning potentials for motion in: b) [0,1] direction 
(a = 0°), c) [-1,1] direction (a = 45°), and d) [-1,2] direction 
(a = 26.6°). Potentials in units of the single CD potential 
well depth (J). 





FIG. 4. Dynamical phase diagrams for drives along [0,1] 
(a = 0°) and [-1,1] (a = 45°). Panels: density-density corre- 
lation functions P(r) at high drives (r = is at the panel's 
center). Nomenclature: MLQ=moving liquid, MSM=moving 
smectic. Others as in Fig. 3. Dotted lines in the pannels for 
B — 2B$ indicate where vortex motion stops. 
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FIG. 7. Dependence of the structure function, 5(k), for 
B = 2B$ and T/T m = 0.83 on driving force, F p , along the 
a — 35° line in Fig. 5. a. Spatial symmetries named as in Fig. 3 
and 4. 




FIG. 8. Dependence of the structure function, S'(k), for 
FIG. 6. Density-density correlation functions, P(r) (r = B= 1.25^ and T/T m = 0.89 on driving force, F P , along the 
is at the panel's center), for low-T infinite- drive phases for a = 40 ° n ne in Fig.5.b. Spatial symmetries named as in Fig.3 
motion in the [0,1] (a = 0°) direction as a function of B. anc [ 4 
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B=2B T/T =0.83 
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FIG. 9. Center of mass velocity direction (0) and com- 
ponent perpendicular to the [-1,1] direction , VI, vs. driving 
force magnitude along lines of constant a in Fig. 5. a. Inset: 
histogram for the distribution of vortices direction of motion 
(0) for a = 35° and F d /F p = 1.5. 
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FIG. 10. Center of mass velocity direction (0) and com- 
ponent perpendicular to the [-1,1] direction ,V1, vs. driving 
force magnitude along lines of constant a in Fig.5.b. 
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FIG. 11. Component of CM velocity perpendicular to 
[-1,1], Vj_, vs. driving force magnitude along lines of constant 
a and at constant temperatures. 
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